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SUMMARY 

Perturbation solutions of equations for laminar incompressible flow- 
in a semiporous channel are presented, and the results are compared with 
those obtained from a fully porous channel. The perturbation parameter 
measures the amount of suction or injection (blowing) at the porous 
wall; positive values denote suction, and negative values denote blowing. 

In the semiporous channel, with a given value of the perturbation 
parameter, blowing decreased the friction parameter at the porous wall 
by 25 percent, and suction increased the wall friction parameter by 50 
percent. For both semiporous and fully porous channels, the pressure in 
the longitudinal direction decreased for zero and blowing values of the 
perturbation parameter. Sufficiently high suction values resulted in a 
pressure rise in the fully porous channel. Either suction or blowing 
has more influence on the local dimensionless velocity profiles and the 
wall friction parameter for the semiporous than for the fully porous 
channel . 


INTRODUCTION 

Experimental results for flow in a rectangular channel with injec- 
tion through a porous wall are given in reference 1. Various porous 
samples were installed in the lower wall; the upper wall was solid. In 
the experiments, conducted with turbulent flow, the injection of air 
through the porous wall markedly influenced the shape of the velocity 
profile. The friction at the porous wall was obtained only qualita- 
tively. The flow was developed with no injection, but not developed 
when there was flow through the porous wall. Although additional tests 
are indicated, it appears that some other geometry may be better suited 
for the experimental investigation, and an analysis is needed to guide 
the selection of the geometry. A significant turbulent- flow analysis 
requires experimental data from the apparatus being sought, but laminar 
flow may be treated theoretically. As part of the study to guide future 
experiments, laminar flow in a semiporous channel is investigated here- 
in. Results are compared with similar results in a fully porous channel. 
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Such an investigation has intrinsic value. Where the flow is re- 
stricted by the proximity of hounding walls, it falls in the category of 
either channel or pipe geometry. Such flow and geometry are found inside 
bodies such as turbine blades and the skins of high-speed missiles whose 
wall temperatures are reduced by the transpiration-cooling process. 
Another application is for the equipment used in the separation of iso- 
topes by the gaseous diffusion process. The present study should be of 
interest in these applications. 

Laminar flow in a channel with fully porous walls was studied ana- 
lytically by Berman (ref. 2). Although only suction through the wall 
was considered, results may be calculated readily for injection through 
the wall. Before the publication of reference 2, Berman pointed out the 
similarity between the flows in porous and semiporous channels (ref. 3). 
The solution, however, was not given for the semiporous channel. 

Berman's results for the fully porous channel, obtained by first-order 
perturbation solutions and, therefore, valid for only smal 1 suction or 
injection, are extended to large values of suction by Sellars (ref. 4). 1 

In the present investigation, a third-order perturbation solution 
for laminar incompressible flow in a semiporous channel is presented. 
Two-dimensional flow is considered. Velocity distributions, wall fric- 
tion, and pressure drop are obtained from the solution for both suction 
and blowing through the porous wall. A second-order perturbation solu- 
tion is also obtained for the fully porous channel for comparison with 
results from the semiporous channel. 


A,B 


g 

_ f _ II 

g > S > 


SYMBOLS 

constants of integration (eqs. (12) and (Al)) 



friction coefficient, p l 
dimensionless stream function (eq. (6)) 

first, second, and third derivatives of g with respect to X 


-*-A different study of the flow in a fully porous channel with large 
suction has been reported recently in an article by S. W. Yuan, entitled 
"Further Investigation of Laminar Flow in Channels with Porous Walls." 
Jour. Appl. Phys., vol. 27, no. 3, Mar. 1956, pp. 267-269. 
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h 

♦ 

channel height 

P 

static pressure 

Apx 

pres sure- change parameter in x-direction (eqs. (27) and (A7)) 

l Apy 

pres sure -change parameter in y-direction (eqs. (28) and (A8)) 

Re 

main-flow Reynolds number, 2u(x)h/ v 

Reo 

main- flow Reynolds number, 2u(0)h/ v 

•^ e w,fp 

wall Reynolds number for fully porous channel, v^h/v 

^ e w, sp 

wall Reynolds number for semiporous channel, v^/2 v 

1 * 

velocity defined by eq. (7) 

•H • U 

fluid velocity parallel to wall 

H 

° V 

fluid velocity normal to wall 

W 

channel width (fig. l) 

X 

distance in main-flow direction 

y 

distance in normal flow direction 

r 

nondimens ional stream function for fully porous channel 

v',r", 

pin 

first, second, and third derivatives of T with respect to A 

A 

nondimensional normal distance for fully porous channel, 2y/h 

X 

nondimens ional normal distance for semiporous channel, y/h 

M* 

viscosity of fluid 

V 

kinematic viscosity of fluid, p/p 

P 

density of fluid 

* T 

shear stress 

¥ 

stream function 
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Subscripts : 
b bottom wall 

fp fully porous channel 

max maximum 

p porous wall 

s solid wall 

sp semiporous channel 

t top wall 

w wall 

x x-direction 

y y-direction 

Superscript : 

average values 


ANALYSIS 

Laminar-Flow Equations 

A sketch of the geometry and flow system for the semiporous channel 
is given in figure 1. The following development parallels that of 
Berman (ref. 2). The channel width is assumed much greater than the 
channel height; therefore, only two-dimensional flow need be considered. 
The Navier-Stokes equations for two-dimensional steady-state incompres- 
sible laminar flow neglecting body forces are (e.g., ref. 5, p. 48) 


du 

du 

1 <^P / 

K2 
a u 

*2 \ 

b u\ 

(1) 

usr- + V"=n — = 
ox oy 

p ox ' 

Vdx 2 by 2 / 

a>v 

dv 

1 dp / 

fb 2 v 

d 2 v\ 

(2) 

^ + ^ = 

" P & +V ( 

^bx 2 by 2 / 


and the continuity equation is 


du bv 

3x + 5y 


0 


( 3 ) 
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The boundary conditions imposed on the system are 
y = 0 (solid wall), u = 0$ v = 0 
y = h (porous wall), u = 0; v = v w = constant 


(4) 


Transformation of Equations 

The transformation of equations (l) and (2) is accomplished by the 
use of the substitutions 




(5) 


and 


g(*) = 


I 


■[“(°, - 


Y_ 

hU 


( 6 ) 


where 


U = u(0) - 


v w* 


(7) 


and u(o) is the average velocity at x = 0. Since g = g(X) only, the 
X will be omitted hereafter. 


The continuity equation (3) is satisfied by the stream function Y, 


since 


and 


1 dl = , 
u = h SX = Ug 


df 

v = " SE = v w« 




( 8 ) 


J 


Use of equations (l) , (2), and (4) to (8) yields 

-!5i = "[r (8g '' - g,2) -S5 8 ''J 


and 


1 dp 

p 3X = 


(v w gg' - h g") 


(9) 


( 10 ) 
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Since the right side of equation (.10) is a function of X only 
(v w ■= constant) , 


S 2 P = 

oX$x 


0 


( 11 ) 


Since U ^ 0, on differentiating with respect to X , using equation (ll)-, 
and integrating, equation (9) becomes 


Re^ sp (g' 2 - gg") + g"' = A 


where 


Re, 




'W,sp - V 

The boundary conditions (eqs. (4)) become 

X = 0 (solid wall), g 1 = 0; g = 0 
X = 1 (porous wall), g 1 =0; g = 1 


( 12 ) 


(13) 


(14) 


Equations (12) to (14), in the present notation, are identical to those 
given in reference 3. 


Solution of Equation 

For small values of Re w ^ gp , the perturbation solution of equation 
(12) is obtained by expanding g and A near Re w sp = 0: 

g = SO + §i^ e w,sp + ®2^ e w,sp + %^ e w,sp + * ' * 

A = Aq + A]_Re w ^ sp + AgRe^gp + AjRe^ gp + . . . (16) 


Substituting equations (15) and (16) into equation (12) and equating 
like powers of Re w ^ sp yield 


S6' ' = A 0 

cr 1 1 ? = A -f. cp rr M — cr 1 2 

S 1 1 s o s o s o 

g s' ' = A 2 - 2g 0 S l + g 0 g l + g l g 0 


(17) 

(18) 
(19) 


g4 


A 3 “ 2g 0 g 2 + g 0 g 2 + g l g l + g 2 g 0 


g 


1 2 


1 


( 20 ) 
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subject to boundary conditions 

for X = 0 (solid wall), = 0; g n = 0 

for X = 1 (porous wall), g^ = 0; g Q = 1, gj^ = 0 

Use of equations (15) to (2l) yields the third-order solution of 
equation (12) : 

g = 3X 2 - 2X 3 + — 6 ^ s -£ (-16X 2 + 27X 3 - 21X 5 + 14 X 6 - 4X 7 ) + 

p e 2 

, (-5,327X 2 - 41, 006 X^ + 206,976X 5 - 243,628X 6 + 

4 : y l y bUU 

99, 792X 7 - 24, 255 X® + 21,560X 9 - 17,248X 10 + 3,136X 1:L ) + 

3 

, 7 ^ (480,343X 2 - 796,742X 3 - 2,659,566X 5 + 6,649,916X 6 - 

_L y ( vDO y I vD4: y \J\J\J 

4,666 , 740X 7 + 2,162, 160X 8 - 2,902,900X 9 + 3,107, 104X 10 - 

2, 255, 526 X 44 + 1,729,455X 12 - 1,234, 800X 13 + 446,880X 14 - 59,584X 15 ) 

(22) 



and 


A = 


- 12 + i Ee w,s P 


2,929 2 

53,900 v,sp 


598,371 3 

147,147,000 w,sp 


(23) 


Absolute convergence of the series given by equations (22) and (23) has 
not been found. The perturbation solution, however, is compared with a 
numerical solution of equation (12) in the section Accuracy of Results. 


* 


Formulas for Velocity, Pressure, and Friction 
The average velocity at an x. location is given by 

f 1 

u(x) = I u(x, X) dX 

Jo 
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Use of equations (6) to (8), (l3), and (l4) shows that 


(24) 


where 


_ Zu(0)h 


Re 0 ~ v 


Then, from equations (8), (24), and (22), 

fu(x,X)l = 6 (x_x 2 ) + - e ^ sp ( -32X + 81X 2 - 105 X 4 + 84X 5 - 28 X 6 ) + 

L U(x) J 70 

2 

ev,sp - ( -5, 327X - 61,509 X 2 + 517,440X 4 - 730,884X 5 + 

2,263,800 v 3 ’ ’ ’ 

349,272X 6 - 97, 020 X 7 + 97,020X 8 - 86,240X 9 + 17,248X 10 ) + 

3 

(960,686X - 2,390, 226X 2 - 13,297, 830X 4 + 

1,765,764,000 v ’ ’ ’ ’ 

39,899,496X 5 - 32,667,180X 6 + 17,297,280X 7 - 26,126,100X 8 + 

31,071, 040 X 9 - 24,810, 786 X 10 + 20, 753, 460 X 11 - 16, 052,400 X 12 + 

6, 256, 320 X 13 - 893,760X 14 ) (25) 

The velocity ratio given in equation (25) is a function only of Re^, sp 
and position X hut not position x. Thus, the flow may he considered 
to he developed. 

Integrating dp(x,x) = ^ dx + ^ dX and using equations (9), (10), 
(12), (13), and (21) result in the pressure 


uA , T v m x 

p(x,X) = p(0,0) + u(0)xjl - 2u( 0 )h 


+ Mh s - — 


(26) 


The pressure change in the x-direction may then he found from equations 
(26) and (23): 





\ 
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Ap s P( x ; X ) - P(°; X ) „ -J*_(_24 + — Re. 
X ' S P i _2/^x hRe 0 \ 


\ pu 2 (0) 


> 929 Re 2 


796^J742 — 3 


35 W ^ S P 26,950 W * S P 
Re, 


147, 147,000 w,sj/ 


V ,SP 

ReQh 


(27) 


Rote that X does not appear on the right side of equation (27) . This 
is a consequence of equation (ll) and results from the assumption of a 
constant v w . Similarly, the pressure change in the y-direction is 

\2 


Ap. 


= p(X,x) - p(0,x) _ p(X,x) - p(0,x)/ 1 


Re, 


0 


y ’ 6P |pvS P pTI 2 (0) \ 2R * 

The shear stress is given by 


w, sp/ 


Re. 


2£l_ . g 2 


w,sp 


(28) 


^ ■ 49 . 

and the friction coefficient may be defined as 

x T 


(29) 


f = 


w 


1 - 2 / , 
2 pu (*) 


(SO) 


Using equations (25), (29), and (30), the friction for the porous wall 
is obtained as 


f sp,p Re - -g"(l) 3 I 19 Re 

8 2 70 w,sp • 4,527,600 sp • 1,765,764,000 ""'w, sp 


+ 124 ' 033 Re 2 m + —2 ^717^1 25 R 3 


(31) 


(The negative sign appears because y increases as the porous wall is 
approached) and, for the solid wall, as 


; sp,s Re _ g"(0) _ ^ 8 

8 2 " 35 W * S P 4 


5,327 p Q 2 480,343 ?a 3 

,527,600 ^sp 1,765,764,000 S P 

( 32 ) 


where Re = 2u(x)h/v. 
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Accuracy of Results 

In order to estimate the accuracy of the solution for the fully 
porous channel, the results from the perturbation solution are compared 
with the results from a numerical solution in table I. The comparison 
is made for Re^sp = -4*0. 

Equation (12) was solved numerically by desk computation using the 
method of Picard (ref. 6, ch. Xl) . The characteristic values necessary 
for numerical integration of equation (12) are g"(0) and A. The initial 
guesses for these values were taken from the perturbation solution. Suc- 
cessive trials finally gave the numerical results listed in table I. As 
is customary in problems of this type, g"(0) and A must be obtained 
more accurately than g' for a given accuracy in g'. 

Comparison of the numerical and perturbation solutions shows that 
g"(0) from the third-order perturbation solution is accurate to within 
0.13 percent, and g"(l) is accurate to within 0.4 percent, (it is 
assumed that the numerical solution is correct.) This difference is 
expected from equations (3l) and (32). A is accurate to within 0.008 
percent. For values of |Re w sp |< 4, the third-order perturbation solu- 
tion should be accurate to within 0.1 percent. 


RESULTS AND DISCUSSION 

The effects of flow through a porous wall are presented and dis- 
cussed for semiporous and fully porous channels. The third-order solu- 
tion for the semiporous channel is given in the ANALYSIS section) the 
second-order solution for a fully porous channel is given in the appendix. 
Some of the pertinent results of these solutions are listed in table II 
for different values of Re Wj S p . 


Velocity Profiles 

The local velocity ratio u(x, X)/il(x) in the channel is shown in 
figure 2(a) for one porous and one solid wall, and for the fully porous 
channel. Both parts of the figure have Re w ^ sp as the parameter. For 
Pe w ^ S p - o = Re w ^fp, both curves are identical and the solutions of equa- 
tions ( 12) and (Al) are exact. 

Blowing (injection) occurs for Re w ^ sp < 0 and is shown qualitatively 
in figure 1) suction occurs for Re w ^ S p>0. Positive values of Re w ^ sp 
place limits on x/ ReQh . These limits are obtained by setting 
n(x) = 0. Then, from equations (24) and (A6a) , since u(0) ^ 0, 
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(33) 


for the semiporous channel, and 


(4k 


max , sp 


1 


(34) 


for the fully porous channel. For |Re w s p | > 0, there is a greater 
effect of Re w ^ S p on the local velocity ratio in the semiporous channel 
than in the fully porous channel. This may he seen hy conqparing both 
parts of figure 2(a) . 

Although u/u(x) is not a function of x, the velocity ratio u/u(0) 
is related to x. From equations (24) and (25), 


These velocity ratios are given in figure 2(b) for Re w S p = 0, 4, and 
-4 with the dimensionless distance x/ReQh as parameter. The growth 
(Re w ^ S p = -4) and retardation (Re w ^ S p = 4) of the velocity ratio u/u(0) 
are evident. These velocity ratios are influenced more by wall Reynolds 
number in the fully porous channel than in the semiporous channel. 

Note that, for a given value of Re w ^ S p, the general shape of the 
curve is unchanged, even though the maximum values are different (fig. 
2(b)). In a porous channel it is evident that fully developed flow de- 
fined by the velocity invariable with flow direction cannot be achieved. 
But, if fully developed flow is defined as a constant value of 
u(x, X)/u(x) for increasing flow direction (with X fixed), then there 
is obtained analytically fully developed flow in both the semiporous and 
fully porous channels when the suction Reynolds number is fixed ( cf . , 
eq. (25) and figs. 2(a) and (b) ) . 


The dimensionless pressure parameter in the y-direction Apy is 
presented in figure 3. These results are obtained by use of equations 
(28) and (A8) . There is no pressure change across the channel when the 



(35) 


and, from equations (A6), 



(36) 


Pressure Distributions 
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walls are solid (v = 0) . In figure 3(a), when there is suction,, the 
pressure parameter Ap y (and hence p(y,x)) is higher than at the solid 
wall over most of the channel. For injection (blowing), the pressure in 
the channel is always less than at the solid wall. For either large suc- 
tion or injection, Ap y ^ sp will always be less than at the solid wall. 
This may be seen by letting Re w ^ gp -*•<*> in equation (28) ) since g 2 is 
always positive, Ap y ^ sp will be negative. 

The pressure-change parameter Ap y ^p p (fig. 3(b)) is referenced to 
the pressure at the center of the channel. Negative values of 
Re w sp = 2Re w pp result in positive values of 4P y ,fp* Again it is 
found that, for either large suction or injection, the pressure parameter 
will be negative over the full channel, except at the center where it is 
zero. Such a solution is given by Sellars (ref. 4). 


The results in figure 3 may be used to estimate the static -pres sure 
drop across the channel. Equations (28) and (A8) may be expressed as 


p(\,x) - p(0,x) L Re w, sp 


| pu 2 (0) 


Re 0 


Ap. 


where Ap y 
that, when 


is given in figure 3. Use of this equation and figure 3 shows 
Ap y is large, Re w ^ sp is small, and vice versa. Hence, for 
the usual value of Re w sp /ReQ (say 0.01), the pressure change across the 

channel cross section should be small compared with the inlet dynamic head 

- 2 / 
pu 


(0) j . This pressure change is also small compared with the local 
dynamic head pu^(x)) when there is blowing. With suction, however, the 
change may not be small compared with the local dynamic head, especially 
for large x/Re 0 h (approaching (x/Re n h) may ) ♦ 


The pressure -change parameter in the x-direction Ap x is given in 
figure 4 for 0 < x/(ReQh) < 0 .25 and in figure 5 for 0 < x/(ReQh) < 0 .035 . 
For both the semiporous and fully porous channels these figures indicate 
that Ap x decreases for Re w ^ sp <0 (blowing). When Re w ^ sp >0, how- 
ever, the analysis indicates a pressure rise Ap x >0 in the fully porous 
channel for sufficiently high values of Re W; sp (a similar result was 
noted in ref. 3). The same situation should obtain if solutions for 
higher values of Re w sp are available for the semiporous channel. An 
explanation for this pressure rise may be as follows: 


When solved for the pressure gradient, a momentum balance for a unit 
width of a channel yields 


Sp _ _ 1 

3x h 



+ T w,b) + 



ay 


(37) 
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The wall shear stresses on the top and bottom walls are positive and of 
the same order of magnitude. When there is fluid injection, the momentum 
integral increases in flow direction; therefore, the pressure gradient 
decreases in flow direction. For the impermeable channel with fully- 
developed flow, the momentum integral does not change in flow direction, 
so that the pressure gradient is still negative. With fluid suction, 
however, the integral term is negative. If the suction is sufficiently 
large, the integral term overcomes the shear terms and a pressure rise 
results (Ap x >0). 


Friction 

The effects on the friction parameter of flow through the walls of 
a channel are presented in figure 6. The friction parameter is given for 
both the semiporous and fully porous channels for -4<Re w ^ S p<4. The 
friction parameter at Re w ^ S p =0 is well established (ref. 7, p. 309; 
ref. 8, p. 51). Since there is a different boundary condition at each 
wall in equation (12), the effects of the flow through the wall are dif- 
ferent at each wall of the semiporous channel. At the solid wall, blowing 
increases the friction parameter; at the porous wall, blowing decreases 
the friction parameter (fig. 6) . The friction parameters at the porous 
wall for the semiporous channel for R e w , S p = 4 and Re^, S p = -4 are 
50 percent above and 25 percent below, respectively, that of an imperme- 
able wall. Corresponding numbers for the fully porous channel are 7 and 
4.4 percent. 

Although the friction parameter fRe/8 is not a function of x 
(e.g., eq. (31)), the friction coefficient f is related to x, since 
Re = 2u(x)h/v. The average velocity u(x) is increased by blowing and 
decreased by suction. The effects of variations in the wall Reynolds 
number Re w S p, therefore, are greater on the friction coefficient f 
than on the friction parameter fRe/8. 


SUMMARY OF RESULTS 

Analytical solutions for incompressible laminar flow in a semiporous 
channel have been presented. A perturbation method was used to obtain 
the solutions. Results were compared with those for a fully porous 
channel. Some of the principal results are as follows: 

1. Either suction or blowing is more influential on the wall fric- 
tion parameter for the semiporous than for the fully porous channel. For 
wall Reynolds numbers of 4 and -4, the friction parameter at the porous 
wall of the semiporous channel is increased by 50 percent and decreased 
by 25 percent, respectively, from the value for impermeable channels. 
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In the semiporous channel, blowing diminishes the friction parameter at 
the porous wall and increases it at the solid wall; suction at the porous 
wall acts conversely. The friction coefficient behaves similarly. 

2. The pressure in the longitudinal direction decreases for zero 
and for blowing values of the wall Reynolds number. This pressure drop 
occurs in both the semiporous and fully porous channels. . A pressure in- 
crease was found for higher suction values in the fully porous channel 
and is indicated for the semiporous channel. 

3. For both the semiporous and fully porous channels, the pressure 
change across the channel is small compared with the inlet dynamic head. 
Compared with the local dynamic head, the change is small for injection 
but may not be small for blowing at large values of the dimensionless 
distance x/ReQh. 

4. For a given value of the wall Reynolds number, the local velocity 
ratios are changed more in the semiporous than the fully porous channel. 
The. velocity ratio built with inlet velocity acts conversely. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 25, 1956 
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APPENDIX - FORMULAS FOR A FULLY POROUS CHANNEL 


Small Suction or Injection 

The equation for a fully porous channel corresponding to equation 
(12) for a semiporous channel is given in reference 2 as, in the present 
notation. 


Re 


V,fp 


(r ’ 2 - it") + r 


= B 


(Al) 


where 


Re. 


v„h Re 


w,fp “ 2v 


= BIT = r = r( A ), A = 2y/ h 


y being measured now from midchannel (fig. 2(a)), since the flow is 
symmetrical about this plane. Boundary conditions are 

A = 0 (midchannel), T = 0; T =0 

A = 1 (porous wall), r = 1; f =0 

The zero- and first-order solutions from reference 2 are 


r 0 =! (3A- A 3 ) 

r = JL (3A 3 - 2A - A 7 ) 
1 280 ' 


B q . -3 

R _ 81 
B 1 “ 35 




(A2) 


and the second-order equation is given as 

nil! 


- b 2 + r 0 q + lys - zr' 0 rl ^ 


subject to 


A = 0 
A = 1 


r 2 " 0 " r 2 


r = 0 = r 

2 2 


J 


(A3) 
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Integration of the equation for T and use of equation (Al) yield 

r = r 0 + r i R %fp + r 2 Ee w,fP = ^ (3A ' A5) + ^ (3A3 ' “ - A?> + 

R 2 

(-4.921A + 6,132A 3 + 1,386 A 7 - 2,695A 9 + 98A 11 ) 

9 , 055 , 200 

(A4) 


2 _ „ . 81 ^ . _ 234 — Re^ fp (A5) 


B - B 0 + B i Re w ,fp + B 2 Re v,fp “ " 3 + 35 Re w,fp ' 137475 


From these relations, the velocity, pressure change, and friction 
may be obtained in the same manner as for the semiporous channel: 


« - (1 - A 2 ) + (9A 2 - 2 - 7A 6 ) 


= r' 


u(x) ■ 2 

Re 2 


280 


'w. 


? f P - (-4,921 + 18, 396 A 2 + 9,702A 6 - 24,255A 8 + 1,078A 10 ) 

9,055,200 v ' 

(A6) 


where 




Ap 


x,fp 


= p(x,X) - p(0,X) = SBvx r- 


\ pu 2 (0) 


[*»■¥] -(fiX; 


u 2 (0)h 2 


1, B72 Re 2 J 1 - 4 R e ^ .m 
13,475 Ke v,fv)\- Re 0 h / 


(A6a) 


oa j_ 648 p 
24 + 35- Re w,fp 


(A7) 


. p(A,x) -p(0,x) p(A,x) -p(0,x)/ Re 0 \ _ 2 fr'-lWO)] 

Ap y,fp “ x dt 2 ± pn 2 (o) \ 4Ee »,fp/ Re »,fp 


2 PV « 


r 

(as) 


f„ Re 


fp -r"m _ « , 3 Rp . 5 j 516 R e 2 ^ 

~8 — “ r ^ “ 3 + 35 Re w,fp + 565,950 R w,fp 


(A9) 
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TABLE I. - COMPARISON OF PERTURBATION AND NUMERICAL 


SOLUTIONS OF EQUATION (12) FOR Re w ^ fp = -4.0. 


Perturbation 


Numerical 
(Picard 1 s method) 

7.756 

g"(o) 

7.74 7 

-21.953 

A 

-21.955 

-4.508 

g"(l) 

-4.525 

0 

8(0) 

0 

1.000 

g(l) 

1.000 

g' 

X 

g' 

0 

0 

0 

0.360 

.05 

.367 

.667 

.10 

.671 

1.125 

.20 

1.129 

1.397 

.30 

1.400 

1.509 

.40 

1.511 

1.486 

.50 

1.488 

1.351 

.60 

1.351 

1.121 

.70 

1.120 

.810 

.80 

.808 

.431 

.90 

.431 

.221 

.95 

.221 

0 

1.00 

0 


TABLE II. - VALUES FOR USE IN FRICTION AND 
PRESSURE- CHANGE RELATIONS 


Wall Reynolds 
number, 

Re w, sp 

Semiporous channel 

Fully porous 
channel 

-g"(i) 

2 

g"(o) 

2 

-A 

-r"(i) 

-B 

4 

4.623 

2.084 

3.786 

3.210 

-1.559 

3 

4.103 

2.311 

5.620 

3.151 

.432 

2 

3.665 

2.540 

7.610 

3.095 

.703 

1 

3.300 

2.771 

9.742 

3.045 

1.847 

0 

3 

3 

12 

3 

3 

-1 

2.755 

3.227 

14.37 

2.960 

4.161 

-2 

2.554 

3.450 

16.82 

2.924 

5.332 

-3 

2.391 

3.668 

19.36 

2.893 

6.510 

-4 

2.254 

3.878 

21.95 

2.868 

7.698 


9 sov 
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Geometry 
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Figure 1. - Geometry and flow in semiporous channel. 
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(a) Local velocity ratio. 


Figure 2. - Velocity distribution in porous channels. 
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Figure 2 


Concluded. Velocity distribution in porous channels. 
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Figure 3. - Pressure change in y-dird'ction for porous channels. 
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(a) Semlporous channel. 


(b) Fully porous channel. 


Figure 5. - Pressure change In main-flow direction for porous channels; 0 < --- * ■ <0.055. 
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Figure 6. - Effects of wall friction on porous channels; 
-4< «e ¥jSp <4. 
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